Electromagnetic fluctuation-induced interactions in randomly charged slabs 
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Randomly charged net-neutral dielectric slabs are shown to interact across a featureless dielec- 
tric continuum with long-range electrostatic forces that scale with the statistical variance of their 
quenched random charge distribution and inversely with the distance between their bounding sur- 
faces. By accounting for the whole spectrum of electromagnetic field fluctuations, we show that this 
long-range disorder-generated interaction extends well into the retarded regime where higher-order 
Matsubara frequencies contribute significantly. This occurs even for highly clean samples with only 
a trace amount of charge disorder and shows that disorder effects can be important down to the 
nano scale. As a result, the previously predicted non-monotonic behavior for the total force between 
dissimilar slabs as a function of their separation distance is substantially modified by higher-order 
contributions, and in almost all cases of interest, we find that the equilibrium inter-surface separation 
is shifted to substantially larger values compared to predictions based solely on the zero-frequency 
component. This suggests that the ensuing non-monotonic interaction is more easily amenable 
to experimental detection. The presence of charge disorder in the intervening dielectric medium 
between the two slabs is shown to lead to an additional force that can be repulsive or attractive 
depending on the system parameters and can, for instance, wash out the non-monotonic behavior 
of the total force when the intervening slab contains a sufficiently large amount of disorder charges. 



I. INTRODUCTION 

Patterned and heterogeneously charged materials, in 
particular if the heterogeneity is disorder induced, have 
received much attention in recent years in a number of 
different research areas. For instance, randomly charged 
polyelectrolytes [l[ and patchy colloids Q show dis- 
tinct collective and thermodynamic properties than or- 
dinary colloids and charged homopolymers 0, Q. Pro- 
teins for instance represents a prime example of biolog- 
ical molecules exhibiting heterogeneous and highly dis- 
ordered charge distributions. The high specificity and 
selectivity of protein-protein interactions is one of the 
fundamental problems of molecular biology and requires 
an understanding of the interaction between randomly 
patterned surfaces @. Another example which has been 
in the focus of recent experimental investigations is the 
problem of interaction between surfactant-coated sur- 
faces which exhibit unusually strong and long-range at- 
tractive forces [(| , shown to stem directly from the pres- 
ence of quenched random domains (patches) of positive 
and negative charges on these surfaces Q. 

In fact, most solid surfaces exhibit heterogeneous 
charge distributions that can be highly disordered as re- 
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vealed by recent Kelvin force microscopy measurements 
Such random charges may result from the surface 
adsorption of charged contaminants and/or impurities, 
while even clean polycrystalline samples display patchy 
surface potentials 0, [l(| ■ The patchiness of the surface 
potential is believed to lead to significantly large effects in 
the experiments aimed at measuring the Casimir-van der 
Waals (vdW) interactions between solid surfaces in vac- 
uum. Indeed, recent ultra-high sensitivity measurements 
have shown the presence of an "anomalously" long-range 
interaction which can easily mask the Casimir-vdW force 
at sufficiently large separations [TTl - fl5j . 

In a scries of theoretical papers [l6l - [20j , the effects of 
quenched monopolar charge disorder in the bulk or sur- 
face of dielectric slabs were investigated. It was shown 
that even a small amount of quenched random charges 
can lead to strong long-range interactions between dielec- 
tric slabs. These interactions were shown to result di- 
rectly from the interplay between the electrostatic inter- 
actions generated by the presence of dielectric discontinu- 
ities (the so-called image charge effects) and the quenched 
statistics of the random charges. It is thus remarkable 
to note that such forces exist even for dielectrics which 
are electroneutral on the average but carry a disordered 
charge component. In this case, net Coulomb forces are 
obviously absent and thus the disordered-induced forces 
directly compete with the Casimir-vdW forces. While 
the latter dominates at small separations, the former be- 
comes substantially large and wins at large separations. 
The previous calculations were however performed only 
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within the classical regime, where only the zero Mat- 
subara frequency contributes to the Casimir-vdW force 
P . HH . Strictly speaking, this approximation would be 
valid above the thermal wavelength (around 7 microns at 
room temperature) [2lj . although the contribution from 
higher-order Matsubara frequencies would in fact dom- 
inate at much smaller separations depending on the di- 
electric properties of the materials (e.g., below about 1 
micron in vacuum and 100 nm in a polar medium such as 
water [22[ ) . The zero- frequency results would be relevant 
for the large-distance regime where the above-mentioned 
anomalous force is observed PP. However, at sub- 
micron separations, it would be necessary to examine the 
quantum effects from the higher-order Matsubara modes 
of the electromagnetic field fluctuations. 

In the present work, we shall thus set out to inves- 
tigate in detail the interaction between two randomly 
charged net-neutral dielectric slabs by accounting for the 
full spectrum of electromagnetic field fluctuations in the 
following two cases: i) the slabs interact across a disorder- 
free dielectric continuum, considering both similar as well 
as dissimilar slabs, and ii) the slabs are separated in gen- 
eral by a dielectric layer which itself may contain random 
quenched charges. In both cases, the results can be com- 
pared directly against those reported previously P, P] ; 
hence, we can determine the effects of the inclusion of 
higher-order Matsubara frequencies, which will be com- 
puted via the Lifshitz formalism P, P] , as well as the 
quenched disorder charges in the intervening medium. 

These results thus generalize the analysis of the dis- 
order effects to all ranges of separation down to the 
nano scale (as long as the continuum dielectric model 
employed within the Lifshitz formalism remains valid). 
We can then draw conclusions regarding the crossover 
between different scaling regimes for the interaction be- 
tween slabs, which were missing from a zero-frequency 
analysis [13, EH- I n particular, we show that the char- 
acteristic ~ D~ x decay P] of the total force with the 
distance D between two randomly charged but otherwise 
(dielectrically) identical semi-infinite slabs sets in well 
within the retarded regime (around, e.g., 50-500 nm). 
Therefore, it is found that the interaction crosses over to 
this ~ D^ 1 disorder- induced behavior from the standard 
(retarded) ~ _D~ 4 Casimir-vdW behavior rather than 
from the classical zero-frequency ~ D~ 3 behavior. It 
turns out that even for highly clean samples (with disor- 
der charge densities down to 10 -9 nm -3 ), the magnitude 
of the disorder-induced force is substantial if compared 
with the Casimir-vdW force. 

For dielectrically dissimilar slabs, we show that the 
non-monotonic behavior P] of the total force as a func- 
tion of distance persists when higher-order Matsubara 
frequencies are included. However in almost all cases, the 
equilibrium separation defined through the zero of total 
force (which can represent either a stable or unstable free 
energy extremum) is shifted to separations that can be 
substantially larger than those predicted within the zero- 
frequency theory (l8| . This is an important consequence 



of our analysis and suggests that the non-monotonic fea- 
tures of the interaction force between dielectric slabs 
could be easily amenable to experimental measurements 
in this regime fl5| . Such non-monotonic interaction pro- 
files have received a lot of attention in the context of the 
Casimir effect in recent years and may arise in the case 
of metamaterials [24[ and/or other exotic materials such 
as topological insulato rs [251 . as well as in certain non- 
trivial geometries P, |27| . In our analysis the ~ D^ 1 
behavior of the interaction force for identical slabs and 
the non-monotonic force profile for dissimilar slabs repre- 
sent characteristic fingerprints of the charge disorder and 
can thus be useful in assessing whether the experimen- 
tally observed interactions in ordinary dielectrics can be 
interpreted in terms of disorder effects. 

The organization of the paper is as follows: In Section 
HT1 we introduce our model and the details of the formal- 
ism employed in our analysis. The results for two semi- 
infinite slabs interacting across vacuum or a disorder-free 
dielectric layer are discussed in Section IIII Al and those 
for the case where the slabs interacting across a dielectric 
layer which itself contains disordered charges is discussed 
in Section lill Bl We conclude our study in Section ITVl 

II. MODEL AND FORMALISM 

We consider a plane-parallel three-slab system consist- 
ing of two semi-infinite regions of dielectric response func- 
tions ei(w) and £2(w) and an intervening slab of thick- 
ness D and dielectric response function e m (ui) (see Fig. 
[TJ. All three slabs are assumed to carry a disordered 
monopolar charge distribution, p(r) [28| . The disordered 
charge distribution is taken to have a zero mean value 
((p(r))) = 0, which ensures that the slabs are net-neutral, 
and a two-point correlation function 

((p(r)p(r')})=g( e -Q';z)6(z-z% (1) 

where ((■••)) denotes the average over all realizations of 
the charge disorder distribution, p(r). Here, g = (x,y) 
denotes the lateral directions in the plane of the slab per- 
pendicular to the z axis where the bounding surfaces are 
taken to be located at z = ±D/2. The above form of the 
correlation function thus implies no spatial correlations 
in z direction and can be thus applicable in general to 
layered materials. In lateral directions, we have a statis- 
tically invariant correlation function whose specific form 
may depend on z, i.e. 

Q{q-q';z)=9{z)c{q-q';z), (2) 

where 

f ci(x) z<-D/2, 
c(x;z) = ^ c m (x) \z\<D/2, (3) 
[ c 2 (x) z > D/2, 

and we shall further assume that the disorder variance 
g{z), which gives the density of random quenched charges 
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FIG. 1: (Color online) We consider two semi-infinite net- 
neutral regions with dielectric response functions ei(a>) and 
£2(^-0 with the static dielectric constants, ei and £2, respec- 
tively, and monopolar charge disorder distributions (shown 
schematically by small light and dark patches) described with 
variances gi and <?2 • The intervening slab may be either vac- 
uum or in general a net-neutral dielectric material of thickness 
D and dielectric response function e m (uj) with the static di- 
electric constant of e m and charge disorder of variance g m - 



in the bulk of the slabs [l7| > ^ s given by 




z < -D/2, 
\z\<D/2, 
z > D/2. 



(4) 



We shall not deal with effects due to disorder in the di- 
electric response of the interacting media [29|, H(| , which 
presents an additional source of disorder meriting further 
study and focus only on quenched disorder (see Refs. 

HH, HH for the cases studied with annealed or 
partially annealed disorder, or with mobile ions on or 
in between the randomly charged surfaces on the zero- 
frequency level). The quenched model is an idealization 
of the real nature of random charge distributions that can 
in general exhibit dynamical behavior, but these effects 
are expected to be small since ion relaxation processes are 
extremely slow on the scale of the Matsubara frequencies. 

We base our analytical calculations on no other as- 
sumption regarding the lateral correlation function, so 
the expressions in what follows may be applied straight- 
forwardly to some rather general cases, such as disorder 
distributions with a "patchy" structure characterized by 
a lateral correlation function decaying over a finite cor- 
relation length (l8j |. Although, for the sake of brevity, 
we restrict our final discussion in this paper to the case 
where the disorder distribution is statistically homoge- 
neous and uncorrelated, c{g — g'; z) = 8(g — g'), thus 



({p(r)p(r')))=g(z)5(r-r'). 



(5) 



In our previous works, we derived the partition func- 
tion of the system defined above for the case where the in- 
tervening medium is free from any kind of charge disorder 



and the electromagnetic field fluctuations are taken into 
account only on the zero-frequency level [13, E3 ■ The 
latter would be a valid approximation only at sufficiently 
large separation distances, D, or sufficiently high tem- 
peratures, T. In the present work, we shall account for 
all higher-order Matsubara modes of the electromagnetic 
field fluctuations, which become increasingly important 
at small separations down to the nano scale. It is easy 
to see that when the disorder is perfectly quenched, as 
indeed we assume here, these charge sources only couple 
to the zero-frequency mode and thus do not mix with 
the higher-order frequency modes of the electromagnetic 
field fluctuations. We do not delve further into the de- 
tails of the derivation of the free energy of the quenched 
system, which can be written, after averaging over vari- 
ous realizations of the charge disorder (see Ref. [l8| for 
details), in an additive form as 



T = TV 



dW 



•7"dis- 



(6) 



The first term on the right hand side above is the 
Casimir-vdW interaction free energy, which is obtained 
in the Lifshitz form of the surface free energy density as 



vdW 



EE ln 



t 



Q n=0 

2D/i m (j£„) 



x e 



+ [(TM) -» (TE)], (7) 



where k-&T =1/(3 and S is the surface area of the slabs. 
The free energy is normalized in such a way that it tends 
to zero at infinite separation distance D and TM and TE 
correspond to transverse magnetic and transverse electric 
modes. In the Lifshitz formula the Q summation is over 
the transverse wave vector and the n summation (where 
the prime indicates that the n = term has a weight of 
1/2) is over the imaginary Matsubara frequencies 



2imk^T 



(8) 



where h is the Planck constant divided by 2ir. All the 
quantities in the bracket depend on Q as well as £ n . We 
have defined 



(TM) 



(9) 



which quantify the dielectric mismatch across the bound- 
ing surfaces between the three different slabs labeled by 
a.,/3 = 1, m, 2. Also n a (i^ n ) for each electromagnetic field 
mode within the slab a is given by 



2 , e a(*£n)Mc(( l £n)£n 



(10) 



where c is the speed of light in vacuo, Q is the magnitude 
of the transverse wave vector, and e a (it; n ) and fJ. a (i^ n ) are 
the dielectric response function and the magnetic perme- 
ability of the corresponding slab at imaginary Matsubara 
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frequencies, respectively. For the sake of simplicity we 
assume that for all slabs /J. a (i£, n ) = 1- Note that e(z£) 
is standardly referred to as the vdW-London dispersion 
transform of the dielectric function and follows as 1331 



2 [ c 
e(%£) = 1 + - / 
n Jo 



to Im[e(u;)] 



duj, 



(11) 



being in general a real, monotonically decaying function 
of the imaginary argument £ [2l|, [22[ . 

For the TE modes everything remains the same except 
that in this case 



«/s(«£n) + «a(«£n) ' 



(12) 



The general contribution from disorder charges follows 
in an exact form as [H[ 

T dis = | fdvdr'g(g- g';z)S(z- z')G(r,r') = 



| J drdr' g(z)c(g- g';z)S{z - z')G(r, r'). 



(13) 



where G(r,r') is the zero-frequency (electrostatic) 
Green's function defined via 



e V • [e(r)VG(r,r')] = -5(r - r' 



(14) 



for the zero-frequency or static dielectric constant profile 
defined as 



e(r) 




z < -D/2, 
\z\<D/2, 
z > D/2. 



(15) 



Equation (|13[) is valid for any arbitrary disorder corre- 
lation function G{q— q'; z) and dielectric constant profile 
e(r). For the particular plane-parallel three-slab model 
considered in this work, the Green's function G(r, r') can 
be calculated from standard methods and one finds 



-h 



S 



dQ 



-2QD 



Ai 



Q l-AiA 2 e 

£l9m ( Til 



-2QD 



(ei +e m ) 2 



(Q) 



(£2 + £ m ) 2 (ei + £m) 2 (£2 + £m) 2 

for arbitrary separation distance D, where 

£?' £m 



Ai 



A,; = 



i = 1,2, 



, (16) 



(17) 



is the static dielectric jump parameter at each of the 
bounding surfaces at z = ±D/2, and 



Pe 2 /(4ire ) 



(18) 



is the Bjerrum length in vacuum (Zb — 56.8 nm at room 
temperature), and c a (Q) is the Fourier transform of the 
correlation function c Q (x). As noted before, we shall fo- 
cus here on the particular case with no spatial correla- 
tions, see Eq. (0, corresponding to c a (Q) = 1 for all 
three slabs a = 1, m, 2. 



Note that Tdis stems from electrostatic interactions be- 
tween randomly distributed disorder charges in the three 
slabs. Due to the dielectric discontinuities across the two 
bounding surfaces, each disorder charge is accompanied 
by an infinite number of electrostatic "images" , which 
are in fact generated by the (static) polarization of the 
slabs. These "image" charges also contribute to the total 
free energy of the system as they interact among them- 
selves and with the actual disorder charges. This type of 
effects are systematically taken into account through the 
electrostatic Green's function G(r, r') and are completely 
included in the above disorder free energy [18 1. 

Our goal is to calculate the effective interaction force /, 
which is mediated between slabs 1 and 2 by both the elec- 
tromagnetic field fluctuations and the disorder charges 
placed in the intervening slab, or equivalently the effec- 
tive interaction free energy T between the two bounding 
surfaces at z = ±D/2, i.e. 

d T 

f = -QD with F ^ D ) = ^dw(D) + F^D), 

(19) 

which can thus be calculated in an explicit form from 
Eqs. and (HH). 



III. RESULTS 

A. Role of higher-order Matsubara frequencies 

In order to bring out the role of higher-order Matsub- 
ara frequencies and compare it with our previous zero- 
frequency results [lU, we shall first proceed by taking 
two slabs interacting across vacuum or a disorder-free di- 
electric medium with g m = 0. We consider three different 
cases: two identical slabs with e\(i£,) = C2(*£) hi vacuum 
and two dissimilar slabs with e\(i£) < e m (i£) < £2(*£) 
and ei(z£), e 2 (»£) < e m (i£). In order to enable a direct 
comparison between these different cases, we assume a 
simple model for the vdW-London dispersion transform 
of the dielectric response functions of the slabs as [22| 



e(zO = 1 



e- 



(20) 



in all cases. This form for the vdW-London dispersion 
transform mimics two characteristic relaxation mecha- 
nisms in the materials (e.g., one due to electronic polar- 
ization and the other due to ionic polarization as is the 
case for Si02 [HI)- The parameters Gi, G2, uii and 0J2 
are chosen such that the required inequality relationships 
between different dielectric response functions at all Mat- 
subara frequencies as well as the desirable values for the 
zero-frequency dielectric constants are obtained. These 
values are given in Table U and are not meant to repre- 
sent any specific material. The qualitative aspects of our 
results do not depend on the particular values for these 
parameters but on the relationships between dielectric 
response functions of the slabs as discussed throughout 
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FIG. 2: (Color online) (a) Magnitude of the rescaled total force, /3\f\l B /S, between two identical semi-infinite and net-neutral 
dielectric slabs in vacuum (e m = 1) is plotted on a log-log scale as a function of the rescaled distance, D/Ib (see Eq. (|19[) h The 
dielectric response function is taken according to Eq. (|20[) with ei = £2 = 3.81, which is appropriate for Si02 (see the text for 
other parameters). The bulk disorder variance is varied in the range gi = 32 = 10~ 6 , 10~ 7 , 10~ 8 , 10~ 9 , 10~ 10 , 10" 11 nm~ 3 (from 
top to bottom). Solid curve shows the pure Casimir-vdW force obtained from Eq. (0). The scaling behavior of the free energy in 
various regimes of separation is shown explicitly, (b) The rescaled total force, fifl% /S between two dissimilar net-neutral slabs 
when the dielectric response functions satisfy the relationship ei(i£) < e m (i£) < £ 2(z£) with the static dielectric constant values 
ei = 5, e m — 10 and ei = 50. Here we assume g m = and 31 = 32 = 10 _6 ,5 x 10~ 7 , 10 -7 , 10~ 8 , 10 -10 nm -3 (from bottom 
to top). Inset shows a closer view of the region around the minimum for gi ~ Q2 = 10~ 6 nm -3 , compared with the results 
obtained by including only the zero- frequency contribution (red curve), (c) Same as (b) but here we take ei(if ), £2(2^) < £m(«C) 
with the static dielectric constant values si = 15, e m = 30 and £2 = 25. 



the text. All calculations that follow are done at room 
temperature T = 300 K. 

Let us first consider the case of two identical slabs 
with equal charge disorder densities, g\ = 92 = ff, in 
vacuum. In this case, we choose the parameters in Eq. 
PU]l appropriate for Si0 2 , i.e., d = 1.098, C 2 = 1.703, 
ui! = 2.033 x 10 16 rad/s, and uj 2 = 1.88 x 10 14 rad/s [H. 
The static dielectric constant of SiC>2 is thus obtained as 
e(0) = 3.81. 



e = e(0) 
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TABLE I: Parameter values for the coefficients Ci and Ci 
in the vdW-London dispersion transform of the dielectric re- 
sponse function (|20[1 are chosen in such a way that the static 
dielectric constants, e = e(0), shown on the left column are 
reproduced. In all cases, we fix the characteristic frequencies 
as ui = 2.033 x 10 16 rad/s and lj 2 = 1.88 x 10 14 rad/s. 



For very low disorder variance, the Casimir-vdW in- 
teraction dominates and thus the total force between the 
slabs in vacuum is expected to follow the standard Lif- 
shitz form for neutral dielectrics and thus go from the 
non-retarded form characterized by the power-law de- 
cay ~ D~ 3 at small separations, through the retarded 
form ~ D~ 4 for larger separations and then back to the 
zero-frequency form which for asymptotically large sep- 
arations scales again as ^ D~ 3 . This behavior is shown 
in Fig [5^. Here we change the disorder variance in the 
range from g\ = g 2 = 10~ n up to 10 -6 nm -3 [Hj], show- 
ing clearly that the disorder effects set in for this range 
of parameters for separations larger than about 50 nm. 
This is well into the retarded regime and is thus beyond 
the regime where the simple zero-frequency results [l8T | 
can be valid. However, once the disorder effects set in, 
they quickly dominate and the interaction force shows 
the characteristic power-law decay ~ D~ x [T3] - Note that 
the magnitude of the total force can increase by orders of 
magnitude as compared with the pure Casimir-vdW force 
(Fig[2k). Also the transitions between various power-law 
regimes may depend crucially on the characteristics of 
the dielectric spectra and may be quite complicated for 
different real materials. 

For identical slabs in vacuum both the Casimir-vdW 
force as well as the disorder-induced force are attrac- 
tive. When the slabs are dissimilar, one may encounter 
more interesting cases where the two effects oppose each 
other [l8| . We now consider a situation where ei(z£) < 
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£ m( l £) < £ 2(*£) and again only the two semi- infinite 
slabs carry disorder charges (g m = 0). In this case, 
the pure Casimir-vdW force is known to be repulsive 
[2TI l36l . [37| but the disorder force may be repulsive or 
attractive depending on the static dielectric constants 
[l8j . In Fig. [2b, we show the results for a case where 
the disorder force is in fact attractive and thus for suffi- 
ciently large disorder variances in the slabs, one obtains a 
non-monotonic behavior for the total force as a function 
of the separation, D. Such a non-monotonic behavior 
is one of the most remarkable features of the interac- 
tions between neutral but randomly charged dielectrics. 
Non-monotonic fluctuation-induced interactions between 
(neutral) materials have received a lot of attention in 
recent years [HHIt]]- The presence of a small amount 
of quenched random charges can thus provide another 
mechanism that can lead to non-monotonic interactions 
when dielectric slabs interact across a dielectric medium. 

Note that there is an equilibrium separation distance, 
Do, where the total force vanishes, which, in the present 
case, represents a stable equilibrium separation between 
the slabs (Fig. [2})). Thus, the disorder effects give 
rise to a bound state between neutral slabs which oth- 
erwise tend to repel each other due to the Casimir-vdW 
forces. On the other hand, we find a maximum attrac- 
tive force at slightly larger separations than Dq. This 
may be used to optimize the thickness of the intervening 
medium in order to achieve the maximum force magni- 
tude between the slabs. The remarkable point is that, 
the zero-frequency calculation [lj| (red dashed line in 
the inset of Fig. [2b) significantly underestimates the 
value of the the bound-state separation; it gives a value of 
Do ~ 1.4Zb — 79.5 nm, while the inclusion of higher-order 
frequencies (black dashed line in the inset of Fig. [2b) 
yields D ~ 3.4i B ~ 193.1 nm for g x = g 2 = 1CT 6 nm" 2 
and the parameter values specified in the figure (inset). 
Hence, a systematic calculation of the Casimir-vdW force 
based on higher-order Matsubara frequencies predicts 
that the above-mentioned non-monotonic behavior can 
occur in a regime which is even more easily accessible to 
experimental verification [l5j than the prediction within 
the zero- frequency calculation [l8| . In the regime that 
the non-monotonic behavior is found, the magnitude of 
the total force in the presence of disorder is typically 
much larger than the pure Casimir-vdW force as may be 
seen from Figs. [5b and c. 

In the case where the dielectric response functions of 
the side slabs are smaller than that of the intervening slab 
ei(«£), £2(2^) < e m (^)i one encounters a situation which 
is opposite to the above case, i.e., the pure Casimir-vdW 
force is attractive [2l[ but the disorder force is always re- 
pulsive and, as expected, dominates at sufficiently large 
separations [l8j |. resulting again in a non-monotonic be- 
havior for the total interaction force between the slabs 
with a maximum repulsive force at sufficiently high dis- 
order variances as seen in Fig. In this case, the sepa- 
ration distance D where the total force vanishes repre- 
sents an unstable equilibrium distance. In other words, 
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FIG. 3: (Color online) (a) The rescaled total force, /3fll/S, 
between two dissimilar net-neutral slabs when the dielectric 
response functions satisfy the relationship ei(«£) < e m (i£) < 
£2(1^) with the static dielectric constant values £l — O, Era — 
10 and £2 = 50. Here we assume g 1 = g 2 = 5x 10~ 8 nm~ 3 and 
change g m = 10~ 6 , 5x 10~ 7 , 10~ 7 , 10~ 8 , 10~ 10 nm" 3 (from top 
to bottom). The red curve shows the results for g m — and 
includes only the zero-frequency contribution in the Lifshitz 
formula. Inset shows a closer view of the region around the 
minimum for this latter case (the point of zero force is shifted 
from Do — 6.31b — 358 nm, red curve, to Do — 9.7/b — 
551 nm, black curve), (b) Same as (a) but for gi — 32 = 0. 



the disorder-induced forces in this case not only oppose 
the Casimir-vdW force but can also give rise to a po- 
tential barrier in the total interaction free energy In- 
terestingly, in this case the location of the equilibrium 
moves to larger values of the spacing as well (e.g., from 
D ~ 1.7l B ~ 96.6 nm to D ~ 2.8^ ^ 159.0 nm for 
gi = 92 = 10 -6 nm -2 and the parameter values specified 
in Fig. [5£, inset). This is directly due to the fact that the 
inclusion of higher-order Matsubara frequencies leads to 
a larger attractive Casimir-vdW force. 
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B. Role of charge disorder in the intervening 
medium 



So far we focused only on cases where the intervening 
slab does not contain any disorder charges and examined 
the role of higher-order Matsubara frequencies. We now 
proceed by examining the effects that may arise from the 
presence of charge disorder in the intervening slab, i.e., 
g m > 0, and compare the results with those we found 
in the previous Sections and with the zero-frequency re- 
sults published elsewhere [l8|, where such effects were 
not included. 

Let us first consider the case where the dielectric re- 
sponse functions fulfill the relationship ei(«£) < e TO («£) < 
62 This situation was analyzed in Fig. when 

there is no quenched random charge in intervening slab, 
i.e., g m = 0. We now increase the value of g m from 10~ 10 
up to 10~ 6 nm~ 3 , while keeping g\ = .92 = 5 x 10 -8 nm 
fixed. As seen in Fig. [3Ji, the non-monotonicity of the in- 
teraction profile fades away as the disorder variance in the 
intervening slab is increased and the interaction between 
the two semi-infinite slabs becomes strongly repulsive. 
This also means that the stable bound-state separation 
Do increases and eventually tends to infinity. It is easy 
to see that this situation arises because adding further 
charges in the intervening slab is energetically unfavor- 
able, although the slab remains charge neutral on the 
average. In order to demonstrate this effect, we set the 
disorder variance in the two semi-infinite slabs equal to 
zero, gi = g 2 = 0. The total force in this case is repulsive 
and increases with g m , Fig. [5b. In fact, one can easily 
show that for the case with g\ = g 2 = 0, the disorder- 
induced force follows from Eqs. © and (fT9|) as 



/3/dis gmlsV 



S 



2e m D' 



where 



'/ 



£]£2 



(ei + e m )(e 2 + e rn ) ' 



(21) 



(22) 



For the case in Fig. [3Jd, we have r\ > 0, and hence the 
disorder force turns out to be repulsive. In general, when 
all slabs carry disorder charges and assuming that gi = 
g 2 = .9, we find 



/3/di, 



where 



X 



9X grnV 

(ei+e 2 ) e m 



, (£2 



£2 + £r, 



£1 



w- < 23 » 



(24) 



The above formulae clearly show that the disorder- 
induced force decays quite weakly as D . 

In the case where the dielectric response functions of 
the bounding slabs are smaller than that of the interven- 
ing slabs ei(«£), £2(1^) < e m(*£)> introducing quenched 
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FIG. 4: (Color online) (a) The rescaled total force, Pfll/S, 
between two dissimilar net-neutral slabs when the dielec- 
tric response functions satisfy the relationship ei(i£), £2(1^) < 
£m(*£) with the static dielectric constant values ei — E2 = 
3.81 (appropriate for SiC>2) and e m = 100. Here we as- 
sume gi — g2 — 5 x 10 
10 _7 ,10~ 7 ,10 

curve shows the results for g m = and includes only the 
zero-frequency contribution in the Lifshitz formula. Inset 
shows a closer view of the region around the minimum for 



nm and change g m = 10 ,5 x 
, 10~ 1U nm~ 3 (from bottom to top). The red 



results (red curve), (b) Same as (a) but for g\ = gi = and 
£1 = 15, e m = 30 and £2 = 25. The red curve shows the 
corresponding zero-frequency results for g m = 10~ 6 nm -3 . 



disorder charges in the intervening slab suppresses the 
maximal repulsive force and shifts the point of zero force 
Do to larger separations as shown in Fig. @Ji and the 
inset. Intuitively, this is because the dielectric images for 
the charges in the intervening slab are of opposite sign 
and thus generate attractions when £i,£2 < £m- Equa- 
tion (|21[) and the results in Fig. HJd for gi = g 2 = 
clearly demonstrates this effect as we have 77 < and 
the only disorder charge contribution comes from those 
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FIG. 5: (Color online) The rescaled total force, j3fl%/S, be- 
tween two net-neutral slabs when the dielectric response func- 
tions satisfy the relationship ei(«f), ea(«£) < with the 
static dielectric constant values ex = £2 = 3.81 (appropri- 
ate for Si02) and e m = 30,40,60 and 100 (from bottom 
to top). Here we assume g\ — — 5 x 10 -8 nm -3 and 
Qm = 10~ 7 nm~ 3 . Inset shows the location of the zero-force 
point, Do, as a function of e m (black dots) compared with the 
data when only the zero-frequency contribution is included in 
the Lifshitz formula (red dots). 



in the intervening slab. 

The above results strongly depend on the dielectric dif- 
ference between the slabs. In fact, the disorder-induced 
force, Eq. ([2"Tj). shows a non-monotonic dependence on 
e m . If we assume e\ = £2 = £ < f m , the attractive 
disorder force will increases in magnitude with e m when 
e m < (1 + v2)e and deceases in magnitude otherwise. 
For instance, in Fig. [5j we choose two identical slabs 
with £1 = £2 = 3.81 but change the static dielectric con- 
stant of the intervening slab from e m = 30 up to 100. 
The total force becomes increasingly more repulsive and 
the distance at which the force becomes zero, £>o, de- 
creases significantly as e m is increased as shown in the 
inset of Fig. [5] (black dots). It should be noted however 
that the zero-force separation is much larger than what 
is expected from a zero- frequency calculation jTsj ] , which 
varies only weakly with e m (inset, red dots). 



IV. CONCLUSION 

We have studied interactions between randomly 
charged but otherwise net-neutral slabs by including first, 
the full spectrum of electromagnetic field fluctuations 
and second, by taking into account the presence of ran- 
dom charges in the medium between the two slabs. This 
is done by employing the Lifshitz theory which includes 
a summation over all Matsubara frequencies, leading to 
a direct generalization of our previous results fl7l Il8| . 
which were obtained based only on the zero-frequency 



effects and were thus valid at large separations. This is 
a crucial new aspect of our analysis of charge disorder as 
it extends the regime of validity of some of the key find- 
ings to the full range of inter-slab separations, especially, 
down to the nano scale (as long as the continuum model 
assumed within the Lifshitz theory remains valid). 

We can thus draw certain important conclusions re- 
garding the role of disorder effects. In particular, it is 
shown that the characteristic ~ D" 1 behavior due to 
quenched disorder [ItJ sets in well within the retarded 
regime (e.g., around 50-500 nm in Fig. ®l) and thus the 
force curves deviate rapidly from the standard (retarded) 
~ D~ 4 Casimir-vdW behavior even for highly clean sam- 
ples with disorder variances down to 10~ 9 nm~ 3 . This 
behavior should be contrasted with the zero-frequency 
theory [lU which predicts that the ~ Z? -1 emerges due 
to a crossover from the classical ~ D~ 3 behavior. 

Another remarkable prediction that follows from our 
present analysis is that the non-monotonic behavior of 
the total force as a function of distance persists when 
higher-order Matsubara frequencies are included; how- 
ever, in almost all cases, the (stable or unstable) equi- 
librium separation is shifted to much larger values than 
predicted based solely on the zero- frequency theory (18} . 
This is important in that it predicts that the disorder- 
generated non-monotonic behavior of the force is more 
easily amenable to experimental measurements that ex- 
pected from zero-frequency calculations. 

We have also shown that the presence of charge dis- 
order in the intervening medium leads to another ad- 
ditive contribution to the total force that can be re- 
pulsive or attractive depending on the system parame- 
ters (such as the relationship between dielectric response 
functions of the media) and can, for instance, wash out 
the non-monotonic behavior of the total force. Such non- 
monotonic behaviors for the interaction between dielec- 
tric slabs have received a lot of attention in the context 
of the Casimir effect and are known to emerge, e.g., in 
the case of metamaterials [24| and/or other exotic ma- 
terials such as topological insulators [25|, as well as in 
certain non-trivial geometries (2(1 [27j . In our analysis 
the ~ £> _1 behavior of the force for identical slabs and 
the non-monotonic force profile for dissimilar slabs repre- 
sent characteristic fingerprints of the charge disorder and 
can thus be useful in assessing whether the experimen- 
tally observed interactions in ordinary dielectrics can be 
interpreted in terms of disorder effects. 

We presented our results explicitly for the case where 
the disorder distribution is statistically homogeneous and 
uncorrelated in space, although the formalism in Section 
HT1 is in general applicable to layered materials exhibiting 
finite lateral correlations for the disorder charges as well. 
Thus, it would be interesting in the future to examine the 
above-mentioned effects in the situation where the disor- 
der distribution consists of random patches (domains) 
of finite size [l8[ , and specifically, when the layering as- 
sumption is relaxed and the correlation domains within 
the slabs have a three-dimensional structure. 
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